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Abstract 

We re-derive the quantum master equation for the decoherence of a massive Brownian particle 
due to collisions with the lighter particles from a thermal environment. Our careful treatment 
avoids the occurrence of squares of Dirac delta functions. It leads to a decoherence rate which is 
smaller by a factor of 2ir compared to previous findings. This result, which is in agreement with 
recent experiments, is confirmed by both a physical analysis of the problem and by a perturbative 
calculation in the weak coupling limit. 
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I. INTRODUCTION 



A classic result of decoherence theory is the rapid decay in the off-diagonal matrix el- 
ements in the coordinate representation of the density operator p(Ri,R 2 ;t) of a massive 
Brownian particle suffering collisions with the lighter particles of a thermal bath. Early cal- 
culations by Joos and Zeh jl| were improved by later authors, and the result of Gallis and 
Fleming Q seems to be the most widely quoted jj. They find, in the limit of an infinitely 
massive Brownian particle, that 

^(R 1; R 2 ;t) 



dt 

where 



-F(R 1 -R 2 )p(R 1 ,R 2 ;t), (f) 



F(K)=en jf ' dqv(q)± J (l - e**-*)*/") |/(gn 2 , qh,) | 2 , (2) 

with m the mass of the bath particles, n their number density, and v(q) dq the fraction of 
particles with momentum magnitude between q and q + dq; hi and n 2 are unit vectors, with 
dhi and dn 2 the elements of solid angle associated with them. The quantity /(q 2 , qi) is the 
scattering amplitude of a bath particle off the Brownian particle from initial momentum qi 
to final momentum q 2 . Gallis and Fleming find e = 2ir. 

We show here that this result is incorrect; the correct result is given by ^ with e = 1. 
Needless to say, this difference does not affect the qualitative conclusion that off-diagonal 
elements decay exceedingly quickly for even macroscopically small | Ri — R 2 1 . Nonetheless, 
experimental techniques are now available that permit the study of the quantum mechan- 
ical loss of coherence by collisions 4]. Therefore not only a qualitative understanding of 
decoherence effects is needed, but a quantitatively correct description is required as well. 
Moreover, derivations of benchmark equations in the theory of decoherence such (1 1121) illus- 
trate the nature of the physics and the assumptions involved, and uncovering the errors of 
earlier results serve as cautionary tales that may facilitate the analysis of situations where 
the physics is more complicated. 

In this paper we present two detailed calculations of the fundamental result (jl|2j) . The 
first is a scattering theory calculation in the spirit of the usual derivations, but one that 
avoids a pitfall of those calculations by using localized and normalized states in the scattering 
calculation. The second is a weak coupling calculation that follows the spirit of master 
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equation derivations undertaken in, e.g., quantum optics. In the first calculation we find (0) 
with e = 1. In the second we find (j2J) with e = 1 and f(qn 2 , gfii) replaced by /B(<?n 2 , <?ni), 
the first Born approximation to that scattering amplitude. This is precisely what would 
be expected, since the second calculation requires the assumption of weak interaction; it 
thus serves to confirm the e = 1 result of the first. Neither of these is the most elegant or 
general calculation one could imagine; the first is rather cumbersome, and the second would 
be neater if generalized to second quantized form Jsj]. But the first has the advantage of 
displaying the physics of decoherence in an almost pictorial way, while allowing a calculation 
involving the full scattering amplitude. And the second, in its simple form, establishes a clear 
connection with the usual approach to decoherence through the master equation approach 
common in quantum optics. Totally separate in their approaches, we feel that together they 
are a convincing demonstration that e — 1. 

These two calculations are presented in sections HP and ITVl below. In section ITTT1 we return 
to the traditional derivation and highlight its inherent shortcomings. We show how it should 
be modified by using a simple physical argument, which leads to a replacement rule for the 
occurring square of a Dirac delta function. This treatment then also yields the result e = 1. 
Our concluding remarks are presented in section 



II. SCATTERING CALCULATION 



To set our notation we begin with a review of the standard approach used to calculate 
collisional decoherence. However, we also wish to point out the difficulties that can arise in 
its application, so we begin in a more detailed way than is normally done. 

To apply scattering theory in a careful way one has to begin with an asymptotic-in state 
|0m) a normalized ket that is the direct product of a Brownian particle ket \<p m ) and a 
bath particle ket The asymptotic-in ket is the result of the evolution of a product ket 
|0m °°' > )|^^ _00 - ) ) at t — — oo to t — under the Hamiltonian that describes the free evolution 
of both particles, without interaction. The effect of the two-particle scattering operator 
S on this asymptotic-in state, S(\4> m )\ijj)), then produces the asymptotic-out state. When 
evolved from t = to t = oo by the non-interacting Hamiltonian, the asymptotic-out state 
yields the actual state at t = oo that evolves from \<pm °°' ) ) |-?/>( _00 )) at t = — oo under the 
influence of the full Hamiltonian. 
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(a) ^ (b) 

FIG. 1: Sketched are the wave packets associated with \4> m ) and \ip) at i = 0. In configuration (a) 
the state \4> m ) \ip) could be taken as both an asymptotic-in state and an initial state at t = 0; for 
configuration (b) that would not be possible. 

In general, of course, \<p m ) does not describe the actual ket at t — that evolves from 
\4>m °°' ) ) |V ;<_00 ' 1 ) at t = — oo, because the evolution of that actual ket involves the particle 
interaction. But if the kets \<p m ) and \ip) are such that the (short-range) interaction between 
the particles has not yet had an effect (e.g., Fig. la but not Fig. lb), then \4> m ) \if)) can be 
taken as the actual ket at t = as well as the asymptotic-in ket. We only consider kets \4> m ) 
and of this form below. 

We now turn to the impending collision of a bath particle characterized by and a 
Brownian particle described by a reduced density operator at t — given by a convex sum 
of projectors |0 m )(0 m |, 

Pin = ^2 Pm 1^ 

m 

= J c£R,i<£R. 2 |Ri) p Q (Ri, R 2 ) (R2I , 

with probabilities p m > 0, J^p m = 1. Here the |Ri,2) label position eigenkets of the 
Brownian particle, and 

Po(Rl ; R2) = ^ Pm (Rl|0m) (0m 

|R*> (3) 

m 

its position representation. Then 

Pr al = Pin®|^)(^| (4) 

can be considered both as the full initial (at t = 0) density operator, and the full asymptotic- 
in density operator. The full asymptotic-out density operator is then 

total o total of 

Pout — °Pin ° 

dR l( lR 2 S (IRx) |V)) Po(Rl, R 2 ) m (R 2 |) Si 



To determine terms such as S (|R) |V'))it is useful to first consider the effect of the S operator 
on direct products |P) |p) of eigenkets |P) of the Brownian particle momentum and eigenkets 
|p) of the bath particle momentum. Since the total momentum commutes with the S 
operator the scattering transformation can be reduced to a one-particle problem, with 

5 (|P) |p)) = y <iq|P - q) |p + q) (-P- ^ P + ^-p-^) , 

where the matrix element here is that of the one-particle scattering operator S Q corre- 
sponding to the two-body interaction acting in the Hilbert space of the bath particle, and 
m* = mMj (m + M) is the reduced mass. In the limit that the Brownian particle is much 
more massive than the bath particle, M ^> m, this reduces to 

S (|P> |p» - J dq |P - q) |p + q) (p + q|<S |p) 

or, moving to a position representation for the Brownian particle, 

S (|R> |p» = J dq |R> e-^' h |p + q) (p + q|<S |p> 

= J dq |R) |p + q) <p + qle-^'^^e'f'^lp) 
= |R> (e- i »" R /' i 5 e i »" R / R |p)) , 
where p is the momentum operator for the bath particle, and so for general states 

S(\R)\i;)) = \R)(e^ h S e^ h \4>)) 

= |R)^ R >, 

where 
and thus 

P ^ = j dR Y dR 2 |Ri)|^ Rl >Po(Ri,R 2 )<^ R2 |(R 2 |. 

Although Po°t al is not the final density operator at t — oo, but only the asymptotic-out density 
operator, it evolves to the final density operator through the non-interacting Hamiltonian, 
and overlaps of the form ^ Ii2 \ip Kl ^ will be preserved during this free evolution. So the final 
reduced density operator for the Brownian particle at t = oo is 

Pfinal = J dR,dR 2 |R 1 )<V R2 |^ Rl >Po(Ri,R 2 )(R 2 | 

= J dR l dR 2 |Ri) p(Ri, R 2 ) (R 2 | , 



where 



p(R 1 ,R 2 ) = <^ R2 |^ Rl >Po(Ri,R 2 ). (5) 

As is well understood, decoherence arises because the bath particle becomes entangled with 
the Brownian particle and the two (asymptotic-out) states |t/> R2 ) and |^ Rl ) resulting from 
scattering interactions associated with the same bath ket and different position eigenkets 
|R 2 ) and |Rx) can have negligible overlap even for |R 2 — Ri| small. The change of the 
Brownian particle's reduced density operator by a single collision is 

Ap(R 1 ,R 2 ) = p(R 1 ,R 2 )- Po (R 1 ,R 2 ) 

= «^|^ Rl )-l)Po(Ri,R 2 ). (6) 

It involves overlap terms of the form 

(^ Ra |^ Rl ) = (^\e- ip - n2,h Sle- ip -^-^ )/h S e ip ^ /h \^) 

= (^SlS^) = tr bath { l S 2 t l S 1 |^)(^|} , (7) 

where the operators 

Sj = e -ip-^^s o e ip ^ /h (8) 

for j = 1,2 are translated scattering operators. We introduce corresponding Tj operators 
according to 

S 3 = l + iT 3 , (9) 
and using the unitarity of the Sj, which follows immediately because S Q is unitary, we find 

SlS ± = 1 + T 2 f 7i - - \t 2 % 



\{r 1 + Tl 



2 ' 2 < 2 



and so 



where 



(V> R2 |v> R = i + (VW), (io) 



4(r 1 + .?)-i(r a + ^). 



Thus the change in the Brownian particle reduced density operator is 

Ap(R 1 , R 2 ) = (VW> Po(Ri, R 2 ) (11) 

The general strategy is to evaluate the matrix element (ip\A\ip) by inserting complete sets 
of momentum eigenstates, 

(V W> = / rfqi^q 2 (VlQa) (q 2 Nqi> (qi|V>> , (12) 

determine (q 2 |.4.|qi), and then perform the momentum eigenstate integrals. Writing S Q = 
1 + i% as well, and using the relations (JBJ) and we find 

(q^lqx) = e 4(qi - Rl - q2 - R2)/ft <q 2 |T V p - (R2 - Rl)/?l T |q 1 ) (13) 
-le^-^-^/^qsl^lqx) 

_l e ,( qi - q2 ).R 2 /^ q2|7 -t To | q ^ 

+ * [^(qi-qaJ-Rl/ft _ e i( qi -q 2 )-R 2 /^] ^ % + T J ^ 

Since the S a operator matrix elements are given by (f| 

(q 2 |5 | qi ) = 5(q 2 - qi ) + ^-5(E 2 - £i)/(cb, qi ), (14) 
where /(q 2 ,qi) is the scattering amplitude, we can identify 

(q 2 |T | qi ) = -L-^-SO/^qi) (15) 

= ^%T /(q2 ' qi) ' 

where i?j = qf/(2m). 

Now in the traditional calculations [110,0] one calculates <9p(Ri, R 2 )/<9t by considering 
the change Ap(R 1; R 2 ) in a time At due to collisions with bath particles that would pass in 
the neighborhood of the Brownian particle, taking the distribution of their velocities from 
the assumed thermal equilibrium of the bath. To calculate Ap(R 1; R 2 ) from one of these 
bath particles, a box- normalized momentum eigenstate, |q) is used in place of a localized 
ket Unlike the |0 m ) states we introduced above, the |0 m ) |q) obviously cannot be 
considered either as asymptotic-in states or as the actual states at t — since the |q) are 
delocalized. Nonetheless, the traditional approach seems to simplify the calculation because, 
as is clear from ()12|) . only diagonal elements (q|^4|q) are required if the limit of an infinite box 



is taken. But from the expression (jl3j) for (q 2 |.A|qi) it is clear that, when a resolution over a 
complete set of momentum states |q') is inserted between Tj and % and the expression (fi3|) 
for the matrix elements of T is used, the diagonal elements (q|*4.|q) involve the square of 
Dirac delta functions S(q — q'). To evaluate these the "magnitude" of 5(0) must be somehow 
set. This is done by relating it to an original normalization volume of the box. While not 
implausible, such a protocol is certainly not rigorous and is open to question. 

To avoid the necessity of this kind of maneuver we will employ bath states that are 
normalized and localized, as is required by a strict application of scattering theory. Before 
addressing the full calculation for a bath in thermal equilibrium we consider scattering 
involving a single state |^). 



A. Scattering of a single bath ket 

From the equations (fT2|) and (fT3|) for (^l^l^) in terms of (q 2 |-A|qi) it is clear that we 
require integrals of the form 

h = J dqidq 2 'u(qi,q2) (q 2 |^ + ^ t |qi) , (16) 
J 2 (R) = J dq 1 dq 2 u(q 1 ,q 2 )(q 2 \Tje l ^ h T \q 1 ), 

which we work out in Appendix ^1 for an arbitrary function w(qi, q 2 ) of the two momentum 
variables. We find that we can write these expressions exactly as 

h = J dqj^dAu(q-^,q+^)M 1 (q,A) (17) 



and 



I 2 (R) = J rfnrfq^dAu(q-y,q+y) e ^ R / ft M2(q,n, A), (18) 

The integration over q covers all momentum space, while A is a two dimensional momentum 
vector ranging over the plane perpendicular to q; n is a unit vector with dh the associated 
solid angle element. Moreover, 

M >«>< A > = 2S5 ('^T-'-T* + na-y-Q+y)) (19) 



and 
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Po ~ ro 




FIG. 2: For this configuration the direct product of p Q and \ip){ip\ can be taken as both the total 
asymptotic-in density operator and the initial density operator at t = 0. 

with 



Q = fi^ 2 + ^- (21) 

With these formulas in hand we can address the expression for (tp\A\ip) once is specified. 
To do this, we take the bath particle wave function (r'|^) to be a Gaussian wave packet 
centered at r„ in position and p D in momentum, 

p ip -(r'-T )/h 

(7ra z ) ' 

and characterized by 

Ax = -=, 
b 

= 71' 

with ab = h. For this minimum uncertainty wave packet we find that the expression for 
(V'lqb) (qi in the integral ((H)) for (tp\A\ip) becomes 

<^|q+ f ><q- y|^> = (_Ly /2 e ^. ro M e -A 2 /(^) e -| q - Po P/^ 

We now assume that this wave packet is located far enough away from the regions of 
space where an initial density operator (jSJ) is concentrated, and with an average momentum 
directed towards the Brownian particle such that the combined density operator (j3J) can be 
taken both as an initial density operator at t = 0, and as the asymptotic-in density operator 
(see Fig. 2). Then using the expressions above we find 

(VW) = (J^j ' J rfnrfqe- |q - p » |2/t2 ^AS(n ) q 1 A) , (22) 



where 



B(n, q, A) = e ^(Ri-R 2 )/VQ-( R 2 - R 1 )/^- iA -( r - R )^ e - A2 /( 4fe2 )M 2 (q, n, A) 

_ 1 -iA.(r -R 2 )/ft e -AV(^) M2(qj fi A) 

+ - [g-iA-Cro-RO/fi _ e -iA.(r -R 2 )/fc] g-AW)^^ A)? 

2 



and where we have put 



R=?±±^. (23) 



The Gaussian functions will keep A within about b of zero and q within about b of p„. We 
now assume that the central momentum p is much greater in magnitude than its variance, 
p ^> b, and hence q 3> b for all q that make a significant contribution; we also assume that 
the scattering amplitude varies little over the momentum range b. Then we can put 

Mx(q,A) » -LRe(/(q,q)) 
M 2 (q,n,A) » -L^ |/(gn, q)| 2 . 

Once these approximation are made the integral over A of the three terms in £>(n, q, A) 
can be done immediately. The integral over A of the first term is not so simple because A 
still appears in Q. In the exponential we have phase factors that vary as 



(24) 



Q (gg ~ Ri) / 2 A 2 n-(R 2 -R 1 ) 

gn-(R 2 -R 1 ) A 2 n ■ (R 2 - RQ 

n 8% 

The first correction term is of order 

6 2 |R 2 -Ri| _ (iRa-R^/a) 

Since q ^> b this term will still be much smaller than unity even if the distance between 
the two positions of the Brownian particle is several widths of the wave packet. We assume 
that |R 2 — Ri| is indeed such this quantity is much less than unity. Then we can replace 
the phase by its leading order expansion 

Q (R 2 - Ri) qh ■ (R 2 - Ri) 

h w h ' 
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in the exponentials of the first two integrals, and the integration over A can be done as well. 
These are two dimensional integrals over a plane perpendicular to q, and so they are of the 
form 

/ dA e- lA -( r °^/ h e- A2 /^ = (27r£) 2 r q (r -R), 

where we have used the fact that ab = H and introduced 

, x exp (- \R 2 - (q • R) 2 1 /a 2 ) , 
T q (R) = lK 1 UJ—L (25) 

which involves R 2 — (q • R) 2 , the square of the component of R that is perpendicular to 
q = q q. In all we find 

(1>\A\il>) = J dd A r °(q), (26) 



where 



A r °(q) = r q (r -R) J dn e l(ci - qiiHRl - R2)/h \f(qn,q)\ 2 (27) 

-i(r q (r -R 1 )+r q (r -R 2 )) / dh |/(gn,q)| 2 
lirift 

+ (r q (r -R 1 )-r q (r -R 2 )) Re (/(q, q)) . 

q 

This is the result we will find most useful when we move to a thermal distribution of bath 
particles. But we close this section with an observation that is of interest in its own right. 
Consider the special case where the size of the bath particle wave packet is much larger than 
the distance between the points Ri and R 2 , 

|Ri-R 2 |/a= |Ri-R 2 |6//l< 1. (28) 

Then in the T q functions of (|27j) we can replace Ri and R 2 by R. Moreover, since the 
integral in (pfi|) restricts q to within a distance of about b of p D , in (|27j) we can replace q by 
p in the scattering amplitudes and in the phase, using the assumption already made that 
they vary little over a range of b; we can also replace the T q functions by corresponding T Po 
functions. The integral in (}2l)j) can then be done, and using (jllj) we find 

Ap(R 1 ,R 2 ) = -p (R 1; R 2 )r Po (r -R) JdA(l- ^(p-^M^-^A) |/( Po n, Po )| 2 . 

The physics here is transparent since 

_ e~ e2 / a2 
r n „(r -R s 



ira 2 
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FIG. 3: A configuration where a> R= |Ri — R2I; I = \/ u 2 — (u • p ) 2 is the impact parameter, 
where u = r G — R. 

where i is the impact parameter of the collision (see Fig. 3). Decoherence occurs only if the 
bath particle does not "miss" the Brownian particle, i.e., if the impact parameter I is smaller 
than the wave packet extension (as limited by the uncertainty principle); its maximum effect 
scales with the integrated square of the amplitude of the normalized bath particle as it passes 
over the pair of points (l/(7ra 2 )); it vanishes as |Ri — R2I — > because the scattering by 
the two points then becomes identical; and it depends only on the scattering amplitude at 
momentum magnitude p a because the variation of that scattering amplitude over the range 
of momentum components included in the wave packet has been neglected. 



B. Convex decompositions of the bath density operator 

To apply the results derived above to a Brownian particle subject to a thermal bath 
we need to describe the effect of the thermal bath in terms of incident, normalized wave 
packets. We start with a single bath particle restricted to a normalization volume Q. In 
thermal equilibrium, at temperature fc^T = the bath state is specified by the density 
operator 

^bath = ^_ e -/3p 2 /(2m) (29) 

provided Q is much larger than the cube of the thermal de Broglie wave length 



V m 

The usual convex decomposition of ()29|) in terms of the delocalized energy eigenstates is 
the obvious one and, aside from the freedom in choosing orthogonal states from among a 
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degenerate set, it is the only one in terms of orthogonal states. But a host of others can also 
found. A particularly convenient set of convex decompositions for our problem at hand can 
be obtained by using 



or, in terms of the pseudo-temperatures k-oT = (3 and k^T = (3 , 

T = f + f . 

In order to use the decomposition (J31)) for p bath ; we write 

e -^(p-p) 2 /(2m) _ e -/3(p-p) 2 /(4m) j fi g-^(p-p) 2 /(4m) 

and take the identity operator in position representation 

In = / dr |r) (r| , 




(31) 



which holds as long as (3 and (3 are both positive and 



1 1 1 



P P $ 



n 



where the integration covers the bath volume fi. We find 




(32) 



where 




3/2 



(33) 



is a normalized momentum distribution function, J £i(p)dp = 1, and the states 




3/4 



(34) 



are characterized by the length scale 
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(compare with (|30p). One then immediately finds 

<rU P > = 2 ^e^'-^e-^ 2 ^ (35) 
so the wave packet |"0rp) is centered at r and has an average momentum p. Indeed, it is of 



the Gaussian form used in the preceding section with minimal uncertainties, b = a/ 2mk B T 

and a = h/b. Thus these wave packets have 

(Ap x ) 2 k B T 
2m ~ 2 ' 

while if we calculate the momentum variance associated with the distribution function /t(p), 



we find 

(S Px ) 2 = k B f 
2m 2 

That is, 

(Affz) 2 (5p x ) 2 _ k B T 
2m 2m 2 

We see that in the class (|32I33|) of convex decompositions of p bath a part of the thermal kinetic 
energy is associated with the size of the wave packets themselves, while the rest resides in 
the motion of the centres of the wave packets. If we take T — > then the wave packets 
are essentially all at rest characterized by a size A — > A, which is the thermal de Broglie 
wavelength. On the other hand, for T ^> T the wave packets are much larger than the 
thermal de Broglie wavelength, and essentially all the thermal kinetic energy is associated 
with the expectation value of the momenta of the wave packets; we have 

(P 2 )- 1 f^,f^„ 3 ; 



2m 2m J ^ ! H 2 ' 
and so a typical speed for the wave packets is 



3k B T 



m 



C. Scattering of a thermal bath of particles 

With the preliminaries of the preceding sections we are now in a position to begin the 
calculation of the effect of a thermal bath on the coherence of a massive Brownian particle. 
We begin with a number of assumptions and choices that will be made, and then discuss 
their applicability and relevance in the context of the calculation. 
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1. Assumptions and choices 



1. We neglect initial correlations, taking the initial full density operator at t = to be a 
direct product of a Brownian particle density operator and a density operator for the 
bath particles in thermal equilibrium, 

p total (t = 0) = Po ®p bath . (36) 

2. We assume that the density of bath particles is much less than A -3 ; then the issue of 
particle degeneracy does not matter and we may consider the density operator of the 
total bath to be just the product of density operators for individual particles. Thus 
we can calculate effects 'particle by particle'. We choose a volume Q much larger than 
any other volume of interest. 

3. We use a convex decomposition of p bath for a single bath particle of the type described 
above, with f < T such that 

f ^ T 

and therefore 

b 2 « (p 2 ) . (37) 

This renders b sufficiently small so that the variation in scattering amplitudes over the 
momentum spread of a wave packet is negligible for essentially all of the wave packets 
in the convex decomposition. 

4. The value of T should also be small enough that the neglect of the variation of the 
scattering amplitudes in the integral is justified, and that we can use the approx- 
imation of neglecting terms on the order of (J24)) above. For the latter we need 

b 2 , 2mk B T 

— R,2 — XM — r R,2 — f*4 ^ 1 

nq nq 

Now for typical wave packets the average momentum p, and hence q, will be of the 
order of mv wp = \J 3mk^T, so this condition becomes 

2tBf 'R 2 -R 1 i^(^ 1/2 ^ |R2 : R ' l «i, 



hv wp V 3 / yTT ^ 

and since T ~ T this reduces to 

T < i— i- ii T 

A 
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5. We choose a coarse-graining time At sufficiently large that 



v wp At > a, 



(38) 



v wp At > |R 2 -Ri|. 



That is, a typical packet travels a distance much greater than its width and much 
greater than the distance between the two decohering sites during the coarse graining 



The second condition is essentially independent of T as long as T m T, and simply 
demands that a typical bath particle wave packet can travel many times the distance 
between the decohering sites during the coarse-graining time. 

It is easy to see that for any Ri and R2 of interest we can meet all these conditions by the 
choice of a large enough At, and we will see below that for a small enough density the change 
in the reduced density operator will be small over any given At. Hence our approximations 
will generally be valid in the low density limit. 

2. Calculation 

The convex decomposition (J32|) we use for the density operator of a bath particle leads us 
to think almost classically about the collision of the wave packets \ip T -p) with our Brownian 
particle. In a naive classical picture the bath particles of a given p that can be considered as 
coming in towards a collision with the Brownian particle in time At are those that lie on any 
of the planes perpendicular to p and extending out a distance pAt/m in the direction — p 
from R (recall (|23|): see Fig. 4). Of course, some of these will completely miss the Brownian 
particle, but none have had a collision with it in the past. For a given p we refer to this 
region of space as R(p). 

Returning to the wave packets, note that those with central positions r close to the Ri 
or R 2 of interest will initially be overlapping with regions of space for which p (Ri, R 2 ) is 
non- vanishing; here any talk of a collision is inappropriate, since at initiation, at t = 0, the 
Brownian and bath particle would immediately be strongly interacting. This is an artifact 



time. Using the expressions for t> wp and a above, and T ps T, the first condition reads 




(39) 
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FIG. 4: For wave packets with momentum p and centers in the hatched region R(p) we take the 
initial density operator to be the direct product p Q ® \4>rp) (V'rpl- For some wave packets, such as 
(A), an actual collision will occur, while for others such as (B) one will not. 

of the assumption (J3HJ) of an initially uncorrelated state, which is clearly unphysical. If the 
uncorrelated state were taken seriously, there would be a rapid "jolt," or shift in the reduced 
density operator due to the set up of correlations [7] . These effects we neglect here, as they 
are implicitly neglected in most such calculations; we do note that, for At large enough, the 
regions R(p) for p of interest will be sufficiently large that only a small fraction of the wave 
packets in the sphere will fall in this problematic class. 

In calculations involving the rest of the wave packets in R(p), we can take the initial 
density operator also to be the asymptotic-in density operator, and use the scattering theory 
calculation for a bath particle wave packet given above. Considering N bath particles in 
volume Q, the total result for Ap(Ri, R 2 ) (recall (fTT| and (}32|) ) is then 

Ap(R 1; R 2 ) = iV Po (R 1; R 2 ) / dp//(p) I ^ (^rp|^|^rp) , 

J J R(p) 

where we assume that the inclusion of the problematic class of wave packets identified above 
will not lead to serious error. Using the result (J2rJ|) from our scattering calculation above, 
we have 

Ap(R 1 ,R 2 ) = np (Ri,R 2 ) / dpp(p) / dr dq A r (q). 

J Jr(p) J (irb 2 y s/ 

If in fact the strong condition ()28|) can be assumed then the integral over q can be evaluated 

immediately, as was done following (J2HJ) . and then also the integral over p can be performed. 
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But this is not necessary. We can simply note that, by virtue of (|3*7j). /t(p) will vary little 
over the range b that e~' q ~ p ' /Speaks and falls. Hence we can replace /t(p) by /t(q) and 
R(p) by R(q), and immediately do the integral over p to yield 



Ap(R 1; R 2 ) =np (R 1; R 2 ) / dq/2(q) / drA r (q). 

J jR(q) 

Since the only r dependence is in the T q , see (|25p. one can now do the r integral for each 
fixed q, putting dr =dr- L dr\ where r" refers to the distance in the direction — q. Since the 
integration over r 1 - is unrestricted in the region R(q) we have 

rfr ± r q (r ± - R^) = 1 

'R(q) 

for Rj = Ri, R2, or R, see (|25j) . On the other hand there is no dependence on r", 



R(q) 



ni 



and so we find 
Ap(Ri,R 2 ] 



np (R 1 ,R 2 ) Jdqfy{q) j dn (e^ h >^- R ^ h - l) \f(qh,q) 



At 

Finally, we recall that T « T and therefore put /t(q) ~ /^(q), where 



c/. equation ()33|). Now if u(q) is the thermal distribution function for the momentum 
magnitude of the bath particles and q = qs, where s is a unit vector and ds the associated 
element of solid angle, we have 

„ (qMq= ^!, (41) 
and hence on a coarse grained time scale we find ()1|2|) with e = 1. 



III. THE TRADITIONAL APPROACH: A REMEDY 

We showed in the preceding section how the problem of evaluating a squared Dirac 
function can be circumvented by expressing the thermal state of the bath particles in an 
over-complete, non-orthogonal basis of Gaussian wave packets (see Eq. ()32|l). However, it is 
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certainly reasonable to explore the possibility of using the standard diagonal representation 
of the thermal bath density operator, which facilitates the formal calculation considerably. 
After all, all the representations of p bath are equally valid and should yield the same master 
equation provided the calculation is done in a correct way It is therefore worthwhile to 
search for a way to deal properly with such an ill-defined object as the "square" of a delta 
distribution function. 

In this section we show how a proper evaluation of the diagonal momentum basis matrix 
elements can be implemented. This leads to an alternate derivation of the master equation 
(ll|2|) . and allows us to highlight the origin of the problem plaguing earlier workers and to 
discuss further implications. However, rather than attempting a mathematically rigorous 
formulation, we base our presentation on a simple physical argument. Our point is that such 
an argument can lead to a prescription for correctly evaluating improper products of Dirac 
delta functions, although this differs from previous naive treatments. 

A. A single collision 

Let us consider again the action of a single scattering event on the Brownian particle 
in position representation p (Ri,R 2 ) and in the limit of a large mass. It follows from the 
discussion in section |H] that after the collision it differs merely by a factor from the initial 
Brownian state, 

p(Ri, R 2 ) = 77(R 1; R. 2 )Po(Ri, R 2 ) (42) 

which is given by 

t?(Ri, R 2 ) = tr bath {e- 4pR2/a 5 t e 4p(R ^ Rl)//i l S e ipRl/ V ath } , (43) 

(see Eqs. (0) and ((7j)). In section ITT1 only pure states p bath = \ip){ip\ of the bath particle were 
considered, but the reasoning is immediately generalized to mixed states. 

The factor r/(Ri,R2) may be called the decoherence function, since it describes the ef- 
fective loss of coherence in the Brownian state which arises from disregarding the scattered 
bath particle. The normalization of p bath implies 

lim 77(Ra,R 2 ) = l (44) 
IR1-R2H0 

which means that the collision does not change the position distribution of the Brownian 
point particle, p(R, R) = p D (R, R). On the other hand, possible quantum correlations 
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between increasingly far separated points will vanish, since a collision may be viewed as a 
position measurement of the Brownian particle by the bath which destroys superpositions 
of distant locations: 

lim r/(R 1; R 2 ) = (45) 

Ri — R.2 |^oo 

This complete loss of coherence implies that the collision took place with a probability of 
one. It could be realized, in particular, by taking the incoming bath particle state to be a 
momentum eigenket in a box centered on one of the scattering sites. 

In thermal equilibrium the density operator (}2T?|) of the bath particle can be written as 

A 3 



p bath = ^^T e -/^/(2 m )| p)(p | 

pePn 
\3 



^-£MP)|P><PI, (46) 



pePn 

with the normalized momentum distribution function ()4()j) at (5 = 1/(&bT). The |p) are 
momentum eigenkets normalized with respect to the bath volume Q, 

and the sum involves those momenta p gPq whose associated wave functions satisfy periodic 
boundary conditions on the box Q. The kets (j4*Tj) form an ortho-normal basis, 

53Sr>5i=2h, (48) 

where 2q is the identity operator in the space of wave functions which are periodic on Q; 
they must be distinguished from the standard momentum kets 

e ipr/h 

<r|p> = (SSp • < 49 > 

which satisfy 

(pIp) = *(p - p') 

and span the full space, 

J dp |p) (p| = J . (50) 
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Since the bath state is diagonal in the momentum representation, an explicit expres- 
sion for the decoherence function (J4"3*|) is readily obtained: 



77 (Ri, R 2 ) -> Jdpfi{p) ^e-^^Sy^-^^So^^'^) 

= j dp(i(p) [l - (p\TjT \p) +e ip < ni -* 2)/h te\Tje ip - {R *~ ni)/h T \p) 
= j dpfi(p) 1 - J dp' (l - ^(p-p'MRx-^)/^ |( p '|T |p)| 2 



In the first line the sum over momenta p GPf7 was replaced by an integral according to the 
usual prescription 



In the second line we introduced the operator T — i(l — S ) and used the unitarity of S , 

K% - To) = , 

as in section |H] and in . The last line follows after inserting a complete set of states 
(|5U|l and noting the relation (|47jl. 

The expression in square brackets in (j5T|) should be well-defined and finite. However, 
it involves two arbitrarily large quantities, the "quantization volume" Q, which stems from 
the normalization of the bath particle, and the squared amplitude of the 7^,-operator with 
respect to (improper) momentum kets. The simple matrix element is given by the expression 

involving the scattering amplitude and a delta-function which ensures the conservation of 
energy during an elastic collision. The squared modulus of (J5*2*j) is not well defined in the 
sense of distributions, but depends on the specific limiting process from which the Dirac 
delta function originates. Yet one would naturally expect that an appropriate replacement 
has the form 

|(p'|T |p)| 2 ^5(p-y)^(p)|/(p',p)| 2 (53) 

with g(p) a function involving the quantization volume. 

We note that (|HT|) displays the correct limiting behavior f)44j) since the expression in 
round brackets vanishes as Ri — > R 2 . On the other hand, for large separations of the points 
the phase in (j5T|) oscillates rapidly and it will not contribute to the integral in the limit 
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|Ri — R2I — > 00. Therefore the limit ()45|) allows to specify the unknown function g(p) in 
(JS3j). One obtains 

fi 1 

</(p) 



(2tt£) 3 a(p)p 2 
with 

cr(p) = j rfn|/(pn,p)| 2 

the total cross section for scattering at momentum p. Formally, this means that one should 
treat the expression involving the squared 5-function and scattering amplitude as 

\S(p-p')f(p'n,p)\ 2 -> 27r ^(p) ^ P ~ P ') l^^^'P)) 2 • ( 54 ) 
B. The master equation 

A master equation can now be derived in the same spirit as above. In the low density limit 
for the bath particles, the decohering effect of each collision can be treated independently 
and the overall decoherence in a short time interval At is determined only by the mean 
number and type of single collisions. 

For bath particles with momentum p the mean number of collisions is given by 
j(p)cr(p)At, where j(jp) = n\p\/m is the flux, and we replace 1/Q by n, the number density 
of bath particles. Hence, the average change in the Brownian state during At reads 

p(Ri,R 2 ) - p (Ri,R 2 ) f, . ,p 



At 



-n Jdp/i(p)^- J dh(l - e *(p-p<0(Ri-R2)/ft) 
|/(pn,p)| 2 po(Ri,R 2 ) 



and noting (JIT]) we have on a coarse-grained time scale which is much larger than the typical 
interval between collisions 

^p(R 1; R 2 ) = -F(Rx - R 2 )p(Ri, R 2 ) • (55) 

with F given by (J2j), again with e = 1. 



C. Interpretation 



It is clear that the derivation of the decoherence function ()43|) does not hold rigorously 
even for volume-normalized ()47|) momentum states, since their amplitude is uniform in space 
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and they cannot be considered as asymptotic-in or asymptotic-out states. Nonetheless, 
the fact that one obtains the "correct" master equation by using the diagonal momentum 
representation (|46|) indicates that it can be reasonable, at least in a formal sense, to extend 
the applicability of (|43|) to volume-normalized momentum eigenstates. 

Then the appearance of the total cross section in the appropriate replacement rule (|54|) 
has a clear physical interpretation. The squared matrix element of the T Q operator with 
respect to two orthogonal proper states may be viewed as the probability for a transition 
between the states due to a collision. The appropriate normalization of the probability 
necessary in the limit of improper states is then effected by the appearance of the total 
cross section a(p) in (JS3J), which is absent in the usual naive treatments of the squared delta 
function. 

This point of view is confirmed by the fact that the rule (p>4"j) . which was derived from a 
simple physical argument (jH3), implies a conservation condition. Integrating ()54|) we have 

(2tt/i) 3 f , ,2 f , . ^ \f{p',p)\ 2 , . 

- J dp |(p |T Q |p)| -> J dp5{p-p) p2o .( p ) = 1 ( 56 ) 

and hence, using (jHOj) and switching to volume-normalized states, 

te|T T tf>->l. (57) 



Inserting the identity (|4"K|) yields 



e i(piwr-i- (ss) 



This is reminiscent of the situation of a multi-junction in mesoscopic physics |8|], or of the 
scattering off a quantum graph where one defines a transition matrix T mn = \t mn \ 2 
which connects a finite number of incoming and outgoing channels. There the t mn are 
the transmission amplitudes between the incoming and outgoing states and the current 
conservation implies 

^ Tmn = 1 With T mn = \t r 



|2 



in analogy to (fHHjl . 

The fact that the conservation relation (fKKJ) has no meaningful equivalence in the contin- 
uum limit VL — > oo is closely connected to the difficulty of evaluating the squared scattering 
amplitude in the momentum representation. It suggests that the diagonal representation of 
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ff can be used in a rigorous formulation of the master equation only if the transition of 
going from a discrete to a continuous set of bath states is delayed until after the square of the 
scattering matrix element is evaluated. A calculation along this line, albeit in a perturbative 
framework, is presented in the following section. 



IV. WEAK-COUPLING CALCULATION 

We now consider an approach that is totally different from the derivation in Section [HI 
Instead of performing a scattering calculation, we obtain a master equation for the reduced 
density operator from a weak coupling approximation that is very similar to the analyses 
of quantum optics. Again the assumption of a low density of bath particles will allow us 
to calculate the effect of the bath particles one particle at a time, so we begin with our 
Brownian particle and a single bath particle restricted to a box of normalization volume fl. 
While we will take the limit Q — > oo in the course of the calculation, we can do it in such a 
way that products of Dirac delta functions never appear. 

In the absence of any interaction between the particles the Hamiltonian reads 

_ V 2 p 2 
H ° = 2M + W 

where m and M are the bath and Brownian particle masses and p and V are their momentum 
operators. The normalized eigenstates of H Q are direct products |P)|p), where |p) is given 
by (JUj) with (J4T?|) . and similarly 

. , iP-R/ft 

(RIP) - -a ■ (59) 



The values of p and P are restricted to a discrete set, p, P G Pn, so that the wave functions 
respect periodic boundary conditions. Our full Hamiltonian is then 

H = H + V(t-K), 

where r and TZ are respectively the bath and Brownian particle position operators, and V 
describes the interaction. 

In the interaction picture the full density operator evolves according to 

ptotai^ = u(t)pf ta \Q)U\t), (60) 
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where U(t) = l+iT(t), and 

T(t) = ~j\ I (t')U(t')dt', 

with 

H^t) = e iH ° t/h V{x-K)e- iH ° t/h . (61) 

Using the unitarity of the time evolution, U\t)U{t) = U(t)W(t) = X, we find from (JtjUj) and 
the definition of T(t) that 

p total (t) _ p total (0) = }-i[ T (t)+T\t),pT a \0)] 

+T (t)pf taX (0)T\t) - ^T\t)T(t)pf tal (0) - ^pf tal (0)T t (t)T(t). 

This equation is exact. The first term on the right side describes a unitary modification 
to the dynamics of the density operator due to the interaction with the bath particles; we 
neglect it here since it will not lead to decoherence. For the other terms we make the 
standard weak-coupling approximations we replace T{t) by 

1 f At 
T = --j^ Hr(t')dt', 

assume an initially uncorrelated density operator, p^ otal (0) = p Q (g> p bath ; (c/. equation (jSEJ)), 
and at time At trace over the bath to find the change in the reduced density operator for 
the Brownian particle as 

A PI = tr bath {T oPo p bath T t - ir o tT oPo p bath - l -p p h ^TlT }. (62) 

To construct T we begin with a Fourier expansion of the interaction potential, writing 

v{x _n) = ^f- ^(q)e iq ' (r ^ )M , (63) 

choosing the prefactor for later convenience. Then forming (|6ip. by inserting resolutions 
of unity in terms of the energy eigenstates between the free evolution terms and the 
potential, we find 

. ^ (2nh? -, Nl rT^ x ^r, fit IP - q| 2 - P 2 

*W = V" £ y(q)|P-q)<P|a«p ( J ^ 

zr |p + q| — p 



X^|p + q>(p|«P ^ )• («) 

pePn V 
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Although this approach leads easily to considering the more general problem of a finite mass 
Brownian particle, we defer that to a later communication. Here we take the infinite mass 
limit for the Brownian particle by taking M — > oo in (|64j) . and we can then write 



P.qePsj 



where 



and 



A Pq =|P-q)(P| 



7q(0 = v(q) E Ip + q)(pl ex P 

Using these in the expression (pj2J) for Ap, we find 



it |p + q| 2 — p 2 
h 2m 



Ap = 2 ^ E (2A P / q /p A P „ q „ - A{,„ q „Ap/ q //9 - p A^„ q „Ap'q') T(q / ,q"), 

P',q',P",q"ePn 



with 



1 /-ZAt /-^t 

T(q', q") = ^ ^ y tr bath { 7q ,(t / )p bath 7 q »(t")} • 
Here we no longer distinguish between the Schrodinger and the interaction picture of the 
reduced density operator since they yield the same evolution in the infinite mass limit. The 
response function Y(q',q") involves the correlator 

'i(t'-t") ( |p + qf-p 2 * 



tr bat 4 7q ^P bath 7 q "(0} = <W \V{4)\ E N ^ 

pGPn 



h \ 2m 



where 

A p = (?|p bath |7) 

is the probability associated with state |p). At this point we include the fact that N particles 
exist in volume Q by multiplying the one particle probability by N. Moreover, for a large 
volume Q we can replace the summation over the bath momenta by an integral, 



pePn 

where n is the bath particle density and p(p) the thermal distribution function ([40)1 . This 
yields 

tr bath { 7q '(^)p bath 7 q »(^)} = S^Q\V(q')\ 2 J np( P )exp (^jp- ^ - ' — 

= (Jq/qZ/fiGq/^-t 7 ), (65) 
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with 

G q (t*-0=n|y(q)| a / Mp)exp 

for all q. Then we have 
(2nhf Q 



i(t' -t") ( |p + q| 2 -p 2 



h 



2m 



dp (66) 



Ap 



a 2 2h 2 



J2 (2A P , q p A P „ q - AJ,„ q A P / q p - PoAp» q Ap. q ) iq(At). (67) 



P',P",qeP r . 



where 



pAt pAt 

Jq(At) = / / G q (t" - t') W . 

Jo Jo 

Now that the Kronecker delta appearing in (|65|) has been summed we can pass to the 
continuum limit Q —>■ oo and switch from the box-normalized states (JoTJj) to the standard 
momentum kets for the Brownian particle, 

e iP-R/n 



(R|P) 



This is done by the replacements 



and 



We obtain 



(2tt£)3 Apq 



(2nhY 



|P-q> (P|=A Pq 



n 



PePn 



dP . 



(68) 



J dP'dP"dq (2A P , q p Ap„ q - Ap, /q Ap, q p - p Ap„ q Ap, q ) I q (At) . 



Similarly, the Fourier transform of (JBUjl yields in the limit 

dr 



via) 



(2Kh) 

(q|V|0) 



:V(r)e 



as the momentum matrix element of the interaction. 
By means of the Fourier expansion 



G q (f) 



duj - 



2tt 



GnMe" 



(69) 



we can write 



7 q (At) = i- y d^G q M 



sin (uAt/2) 
uj2 
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where 



- 2immh ,- ,2 / (3 V' 2 
G n {uj) = \V{q}\ ( 77— ) exp 



Q 



2nm 



f3mh 2 
2q 2 

as we show in Appendix [BJ It follows that for time intervals 



At > (3h 



2 \ 2' 



CO 



2mh) 



(70) 



(71) 



we can set G q (u) ~ <5 q (0), and / q (At) ps G q (0)At for essentially all q of importance in f|UTj) : 
a coarse graining time long enough that (J71|) is satisfied is sufficient to guarantee that Fermi's 
Golden Rule holds and energy conservation is satisfied for the scattering bath particles. We 
then obtain a master equation for the reduced density operator of the Brownian particle, 



^ = 47i 3 h J dP'dP"dq (2A P , q pA P „ q - A P „ q A P , q p - pAp„ q A P , q ) <5 q (0). 

Although not explicitly in Lindblad form this equation yields a completely positive evolution 
of p and can be put in Lindblad form, since G q (0) > 0. 

To see the physics of this result and to compare it with earlier work we go into the 
coordinate representation, p(R 1; R 2 ;t) = (Ri|p(£)|R 2 ), were we find 

dp(Ri, R 2 ; t) 



dt 



— F(Ri — R 2 )p(Ri, R 2 ; t). 



with 



F(R) = 8n 3 h J dq (l - e^'*^) G q (0). (72) 

While this in practice might be a useful expression for evaluating -F(R) directly, to compare 
with earlier work it is easiest to return first to and formally determine G q (0). Inserting 
in (J721) leads to 



F(R) = IQix^tfn dq dp ^t(p) \V(q)\ 2 [l - e" iqR/?l ] 6 



1 1 1 2 2 
|p + q| - p z 



2m 



and with the introduction of a new variable p' = p + q to replace q, 




F(R) = 16vr 4 n 2 n / / rfp'rfpp(p) \V{p' - p) 



1-e 



i(p-p')-R/fr 



(pf-p 2 

2m 



5 



The p' integration can be done in polar coordinates, p' = p'n, i.e, dp' = (p 1 ) dp'dh. Then 
since 

^^H&'->)=>->. 
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we have 

F(R) = 16n 4 h 2 nm J J dndpfi(p)p \V{ph — p)| 2 [l — e *(p-P*W«] , 

which, with the aid of (jUJ, yields 

F(R) = A^tfmn J J pv{p)dpduds \V(ph-ps)\ 2 [l - e *G*-i*)-H/ft] . 

Finally, we use the fact that the first order Born approximation for the scattering amplitude 
is given by [f| 

/B(p',p) = -(27r) 2 m^(p'-p) (73) 

to find 

F(R) =nj J ?-v(p)dp^ \f B (pn,ps)\ 2 [l - e^'^ n ] , 

in agreement with (J2J) for e = 1, if the full scattering amplitude is replaced by its first Born 
approximation. 

V. CONCLUSIONS 

In this article we presented two detailed derivations of the quantum master equation for 
a massive Brownian particle subject to collisions with particles from a thermalized environ- 
ment, and a third argument for the master equation based on physical motivations. They 
represent rather different strategies of dealing with the principal problem that arises in the 
formulation of collisional decoherence: The momentum eigenstates that are the most natural 
for describing the density operator of the thermal bath do not constitute proper states for 
the application of scattering theory. 

By representing the thermal bath with an appropriate over-complete basis, we obtained 
in section |H] the full master equation in a calculation that is mathematically and physi- 
cally convincing, but cumbersome. The perturbative treatment in section II VI permitted the 
avoidance of improper states up to a point where they could be handled in an unambiguous 
way. Yet this calculation hides some of the physics involved, and does not in itself suggest an 
immediate generalization to higher orders in the perturbation parameter. More straightfor- 
ward, but also most ventured, is the use of the replacement rule put forward in section lllll 
It yields the master equation immediately, but lacks the mathematical foundation one would 
wish for an equation describing a process as fundamental as decoherence by collisions. The 
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strategies developed here lead to natural generalizations, both to the more general problem 
of a Brownian particle of finite mass, and to other collisional decoherence processes. We 
plan to turn to them in future communications. 

Different as the three approaches may be, they all suggest that previous results in the 
literature predict decoherence rates that are quantitatively too large. This conclusion is not 
only of theoretical interest. A recent experiment j^J on the decoherence of fullerene matter 
waves by collisions with background gas atoms was sensitive to the presence or absence of 
the factor of 2tt that has been the focus of this article. And indeed, the observed decoherence 
rates are in full agreement with the equation derived in this article, and exclude the previous 
results. 
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APPENDIX A 

In this Appendix we confirm the results (jl7H8|) . beginning from the definitions (|16p. We 
begin by inserting a complete set of momentum eigenstates to write the matrix element 
appearing in I 2 (R) as 

(^teW^lqi) = |rf q '< q2 |T teW|q')(q'|T | qi ) 

= J rfq'e^'- R /' l <q 2 |T t| q '>(q'|r |q 1 ). 

Then using the expression (J15|) for the matrix elements of T and, noting that we can write 
8{q 2 - q')S{q' - q x ) = 5(q 2 - qi)6(q' - q{), we find 

h = J ^0«(qi» <b) (/(*, qi) + /* (qi> q*)) *(* - a), (ai) 

/ 2 (R) = f d A q tT dq ' u ^ qOe^VV, q 2 )/(q / , qi W - gi)%2 - gi), 

and so S(q 2 — qi) appears in both these expressions. They are simplified if we make a change 
of variables from (cji,^) to (q, s), where <\ 2 = q +s/2 and c[i = q — s/2. The Jacobian of 
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this transformation is unity, so <iqi<iq2 = dqds, and 



()(«I2 - <n ) d [ \jq 2 + — + s- q- yg 2 + — - s-qj. 

We do the s integral first, and we single out the component of s along q by writing 

s =z/q + A, 

where A is a two-dimensional vector lying in the plane perpendicular to q. Then ds = dud A., 
and 

<%2 - qi) = S(g(u)) 

for fixed A, where 



, A 2 + u 2 A 2 + v 2 

g[y) = \lq 2 + + qv- \ q 2 + 



4 

which as a function of v vanishes at v = 0. Thus 

%M) = — ^ — = 

\dg/dv\v=o 
and for any function t/(qi, cfa) we have 



4 



qv 



S 2 + ^ 



dq 1 dq 2 U(q 1 ,q 2 )S(q 2 - <7i) 
dqdudA U(q - |, q+ |)<%(v)) 

<iqdA^f/(q-^,q+^), 
g 2 2 



(A2) 



where we have put 



Q 



A 2 


A 




A 


+ T = 


q 

M 2 




q+ y 



qi = q2 



cf. equation ()21|). Using this in the first of ()A1|) we immediately find (fTTj) ; using it in the 
second of flM]) we find JED 

= / SI i"4*»4«4> W(Q,c4)/(Q, q 4), 

where to get from the first to the second line we have put dq' = (q') 2 dq'dh, where dh is an 
element of solid angle, and defined Q as in (|2*T|) . 
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APPENDIX B 



Here we confirm the result (|70p. From the inverse transform of ()69|) we have, using (|66|). 

G q (w) = y dte iW 'G q (t) (Bl) 



= 7r/j, |V^(q)| 2 n y dpdaS (s u (a)) 

where 

2pqa + q 2 

Su{a) = fio;. 

2m 



To get from the second to the third line of (jBl|) we have put 

|p + q| 2 -p 2 2p-q + g 2 

2m 2m 

2pqa + q 2 

2m 

where a = cos 8, and 9 is the angle between p and q, and used (JHJ, writing ds = dad<ft, 
where <fi is the azimuthal angle around q. The function s w (a) has a single root 

2mhu — q 2 
2pq 

which, for there to be a contribution to the integral in ([Bip . must satisfy — 1 < a Q < 1. 
Since p and g are both positive this implies 

\2mhuj — q 2 \ 

P > Pent = 7; • 

2q 

Finally, noting that 

m5(a-a ) 

dis^a)) = , 

pq 

we obtain the expression 

- ixmh , - |2 f°° u(p) 
G q (cu) = V^(q) n dp. 

In thermal equilibrium, 

(3 N 3 ' 2 



v ( p ) = 4tt ( I e -^/(2m) 2 

1 27rm ' 
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this integral immediately yields (|70|). 
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